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Abstract 

Let E{B, Z,N) denote the ground state energy of an atom with 
electrons and nuclear charge Z in a homogeneous magnetic field B. We 
study the asymptotics of E(B, Z, N) as i? — > oo with and Z fixed but 
arbitrary. It is shown that the leading term has the form (In B)'^ e{Z , N) , 
where e{Z, N) is the ground state energy of a system of A^ bosons with 
delta interactions in one dimension. This extends and refines previously 
known results for A^ = 1 on the one hand, and N, Z ^ oo with B/Z^ — > oo 
on the other hand. 



1 Introduction 



The effects of extremely strong magnetic fields (order of 10^ Gauss and higher) 
on atoms and molecules are of considerable astrophysical as well as mathematical 
interest and are far from being completely understood in spite of many theoretical 
studies since the early seventies. We refer to | LS Ya | and | RWHG ] for a general 
discussion of this subject and extensive lists of references. An atom (ion) with A^ 
electrons and nuclear charge Z in a homogeneous magnetic field B = (0, 0, B) is 
(in appropriate units) usually modeled by the nonrelativistic many-body Hamil- 
tonian 



N 



H 



B,Z,N 



N 



i<j 



(1) 
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Here Xi G are the postitions of the electrons, i = 1, . . . , N, A{x) = |B x x is 
the vector potential, and 

/fA = [(iV + A(x))-o-]2 (2) 

with cr the vector of Pauli spin matrices. The Hamiltonian Hb,z,n operates on 
the Hilbert space TCn = /\^L^(R'^; C^) appropriate for Fermions of spin 1/2. In 
this paper we are concerned with the ground state energy 

E{B,Z,N) = ini spec Hb,z,n 

= M{{^,Hb,z,n^) : ^eCo°°(M=^^;C2'')n7^iv,||^||2 = l},(3) 

more precisely the B oo asymptotics of this quantity. Such an asymptotic 
study is relevant at the field strengths prevailing on white dwarfs and neutron 
stars. 

Previous investigations of the asymptotics of E{B, Z, N) have either dealt 
with the case = 1, i.e., hydrogen-like atoms [|AHS|| , or the case when Z and 

tend to oo together with B | [LS Ya] |, | [LSYb| , |^. The most complete rigor- 
ous treatment of the = 1 case so far is ||AHS|| where the following B —>■ oo 
asymptotics was derived: 

E{B,Z,1)/Z'= - i[ln(i?/2)]2 + [ln(i?/2)lnln(i?/2)] 
- [{C + In 2) ln(5/2)] - [In HB/2)]^ 
+ 2(C- l + ln2)lnln(S/2) + 0(1), (4) 

with a constant C (Euler's constant/2). The basic results on the N, Z ^ oo case 
were obtained in [ [LS Ya| ] and | [LS Yb|| . In particular, in ||LS Ya|| it was shown that 
if iV, Z ^ oo with A = N/Z fixed, and B/Z^ oo, then 

f -i^ + l^'-^^' if ^<2 
E{B,Z,N)/{Z-'[\n{B/Z')r)^{ (5) 

I -1 if A > 2. 

The fact that the right side of @ decreases with increasing N/Z as long as 
N/Z < 2 shows that in the limit Z ^ oo, B/Z^ ^ oo an atom can bind at least 
2Z electrons. In [|| some higher order corrections to the leading asymptotics for 
the energy are discussed. 

The main result of the present paper is a derivation of the leading term in the 
B ^ oo asymptotics of E{B, Z, N) where Z and are fixed, but arbitrary. The 
precise statement is as follows: 

1.1. THEOREM (High field limit of the energy). For each fixed Z and N 



(6) 
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where e{Z, N) is the ground state energy of the Hamiltonian 

N N 

hz,N = J2 HVdz^ - Z6iz,)) + J2 ^(^^ - (7) 

i=l i<j 

of N bosons with 5-interaction in one dimension, defined in the sense of quadratic 
forms as 

e(Z, N) = inf{(^, hz,N'^) : ^ G Co°°(M^), W^h = !}■ (8) 

It is trivial to compute e{Z, 1) = — Thus (^) generalizes the first term 
in the expansion (^) to the case > 1. The relevance of the 5-function model 
for the ground state of hydrogen in strong magnetic fields was noted already in 

i- 

We also verify that the mean field limit of e{Z, N) agrees with (|^): 

1.2. THEOREM (Mean field limit). IfZ,N ^ oo with X = N/Z fixed, then 

-iA + |A2-iA3 if A<2 
e{Z,N)/Z-'^{ (9) 
-i if A > 2 



Taken together, Theorems and |1.2| lead to the same high B, high Z limit 



as Theorem 1.4 in ||LSYa|| , where Z ^ oo and B/Z^ oo simultaneously (the 
"hyper-strong" limit.) 

We now describe briefly the strategy for the proof of these results and intro- 
duce some notation that will be used throughout. The first step in the proof 
of Theorem |L1] is a reduction to the subspace H% C Hn generated by wave 



functions in the lowest Landau band. Let 11^ denote the projector on TC%. (Its 
integral kernel is given by Eqs. (p^-(^) in Section 5. Note that 11^ and TC% de- 
pend on B.) Let i?™°^(_B, Z, N) denote the ground state energy of II%Hb,z,n^n- 
It is clear that 

E{B,Z,N) < E^°''\B,Z,N), (10) 

and by Theorem 1.2 in 



E™'^*(5, Z, N) < E{B, Z, N){1 - 5{B, Z, N)) (11) 

where 6{B, Z, N) ^ for i? ^ oo with Z, N fixed. Hence it suffices to prove (P) 
with E{B,Z,N) replaced by E'''"'^ (B , Z , N) . We note in passing that ([Tl|) also 



holds for bosons. In fact, it will become evident in the sequel that Theorem |LT 
is independent of the statistics of the particles. 

To study E'^°^^ {B , Z , N) the next step is to introduce a Hamiltonian for the 
motion parallel to the magnetic field with the coordinates perpendicular to the 
magnetic field as parameters. We write the variables Xj G 

where xj- G and 2;j G M are respectively the components perpendicular and 
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parallel to the field. Moreover, we write {xi, . . . ,xn) = with x"*" = 

(xi, . . . , x]^) e R2A^ and ^ = (^1, . . . , ^^v) e M^. 

In the lowest Landau band the part of (Q) associated with the motion per- 
pendicular to the field is exactly canceled by the spin contribution and only the 
part corresponding to the motion along the field remains. Hence 

^%Hb,z,n^% = ^%Hz,N^% (12) 

with 



N / ^ ^ TV 

^^^^ 



The operator (|T3D contains no derivatives perpendicular to the field and hence 
the variables x-^ can be regarded as parameters for a differential operator in the 
variables parallel to the field. For each x-*- such that all different 

from zero, we consider the one-dimensional Hamiltonian 



1=1 



N 

■2(ja\ _ T2(m>N\ 



acting on (^L (M) = L (M ). The expectation values of Hz^n can be written as 

(VI/, Hz^n"^) = j {^{x^, ■), HzA^^)^^^^^ ■))l2(r-) dx^. (15) 

The next step is a scaling of the variables. In the lowest Landau level the 
characteristic length in the directions perpendicular to the field is B'^^"^. One 
can therefore expect that for the computation of E^°'^^{B, Z, A^), i.e., the infimum 
of (1^) over (normalized) \1/ € H^,, the properties of Hz,n{x_^) for |x^| ~ 5"^/^ 
are decisive. Anticipating this, it is natural to make a transformation of variables, 
{x^-i^ ~^ (B^^^Hl^^ L{B)z) where the scale factor L(B) in the direction of the 
field has still to be specified. The corresponding unitary operator on L^(]R^) is 

U^{z) = L{By/^^{L{B)z), (16) 

and the Hamiltonian transforms in the following way: 

U-^Hz,n{x^)U = L{Bfhl^{B'/'x^) (17) 

where 

N 

hzAt) = E (-^' - ^^B,lyA^^)) + E ^B,ly--yt\i^^ ' (18) 
1=1 i<j 



BSY 19/11/99 



5 



and the potential VB,r{z) is (for r > 0) defined as 

VbA^) = L{B)-\B-'L{Bfr' + z^)-^!\ (19) 

Let -E'z,Af(ai"'") and ef ^(y-*-) denote the ground state energies of Hz,n{x'^) and 
nIv^) respectively. In order to avoid discussions about the domains of the 

Hamiltonians, which in fact depend on whether some of the parameters x:^ (resp. 

y:^) coincide, we define the ground state energies in terms of quadratic forms in 

the same way as (|): 

EzAsi-^) = inf{(^,i/z,iv(x^)^) : ^eCo°°(M^), 11^112=1}, (20) 
el^iy-') = inf{(vI/,/.|^(y^)xl/) : v]/ e Co°°(M^), H^lh = l}- (21) 

These energies are connected by the scaling relation 

EzAB-'V)/L{Bf = el^{y^). (22) 

In the next section we show that with the choice L{B) ~ Ini? the potential 
Vb^z) converges for each r > in the sense of distributions to the delta function 
as -B ^ 0. This is the heuristic basis of Theorem |0|. Since the convergence is 
not uniform in r, however, more is needed for a rigorous proof. In particular, one 
needs estimates on the r-dependence of the convergence VB,r{,z) — ^ ^{z). These 
estimates, stated in Lemmas p.l| and p.2| in the next section, can be regarded as 



variants of Propositions 3.3 and 3.4 in ||LSYa|] and the Appendix in PY| , adapted 



to the problem at hand. They are included here for completeness. 

The upper bound on the energy, given in Section 3, is a straight-forward 
variational calculation. The lower bound is more subtle. An important ingredient 
needed is the superharmonicity of the energy Ez,n{x'^) in the variables xf. This 
result, established in Theorem O, generalizes a corresponding result (Proposition 
2.3) in |P^SYa|| . Superharmonicity implies that the lowest value of Ez^N{B~^^'^y-^) 



for \y^\ > e with e > is obtained at the boundary of the variable range, i.e., 
when either \y:^\ = e or \y:^\ — >■ oo. Variables tending to infinity can be ignored, 
since VB^r{z) — ^ for r oo, so by this result one may in (|T5|) restrict the 
attention to wave functions localized where < (const. )i?^^/^. On the other 
hand, the requirement that only wave functions in the lowest Landau band are 
taken into account in ([T5|) plays the role of a 'hard core condition' that prevents 
collapse, since such wave functions cannot be concentrated on shorter scales than 
0(5"^/^). This statement is made precise in Lemma [57^ . 



The lower bound is obtained in Section 5 by combining Theorem [4.3| , Lemma 
|57^ and the convergence of the potentials VB,r discussed in Section 2. It is 
noteworthy that this lower bound holds also for bosonic statistics while the up- 
per bound holds for fermionic statistics, so that altogether the convergence of 
E{Z, N, B)/(\iaBy to e{Z, N) is independent of the statistics. 
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In section 6 we discuss the delta-function model (0) and in particular prove 
Theorem |1.2| . In the course of the proof we compare (|^) with another model, 



whose ground state energy can be explicitly calculated. This model provides an 
upper bound for the ground state energy of (|^) and has the same mean field limit. 
The Hamiltonian for this model is 



^ 1 
hz,N = J2ip"- '^(^^)) + ^ E '^(l^^l - I^^D- (23) 



i<j 



An interesting feature of this model is the fact that the maximal number A^^c 
of electrons that a nucleus of charge Z can bind is exactly the largest integer 
satisfying 

N^<2Z + 1. (24) 

(This fact is unrelated to Lieb's upper bound for the maximal negative ioniza- 
tion of atoms that does not apply to the Pauli Hamiltonian with a homogeneous 
magnetic field.) A corresponding statement for the Hamiltonian (|^) is not known, 
except in the mean field limit, cf. Theorem 1.2. In this connection it should be 
mentioned that an estimate of the form form Ac < 2Z + 1 + (const.) B^^'^ has 
been derived in [Bn| for a Hamiltonian of a similar type as (|T^). 



2 The high B hmit of the Coulomb interaction 

We define the scaling factor L{B) in the potential (ITW) as the solution of the 
equation 

B'/^ = L{B)smh[L{B)/2]. (25) 
Since J^ia'^ + z'^)'^/'^dz = sinh^"'^(l/a), we have with this choice 

/ VB,r{,z)dz = 1. (26) 

J |2|<r 

for all B. Note also that 

L{B) = In B + 0{\n\nB) (27) 

as -B — s> oo. 

Let ip ^ H^(^) = {"^ • / / Idip/dz]"^ < oo}. Every such ip is a. continuous 
function on M. 

2.1. LEMMA (Delta approximation, part 1). 

1^(0)1' - j VBAzMiz)\^dz\ < L{B)-^ [Ar-i + SX^'^T^'/^r^/^] (28) 
with\ = T = J\d^/dz\\ 
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Proof. It suffices to take r = 1, for the general case follows by scaling z — ^ rz. 
Write the difference on the left side of as Ai + A2 with 



A, 
A2 



VBMm^)\^dz, 



z\>l 



z <1 



The missing term 



A, 



VB,i{z)dz 



z\<l 



\m\' 



is zero because of (^61) . Since |Vb,i(2;)| < L{B) ^ for 1^1 > 1, we have 

< XL{B)-\ 

For \z\ < 1 we have in any case 

\Vb,i{z)\<L{B)-'\z\''- 

Moreover, 

' Kz)\'-\m\' < i^(^)-v^(o)i [1^(^)1 + 1^(0)1] 



(29) 
(30) 

(31) 

(32) 
(33) 



< 



■ 2 



dz' 



< \z\^/^T^/^2\^/^T^/^ = 2A1/^T3/4|; 



1 1/2 



(34) 



Hence 



< 2LiB)-' ( [ \z\-'/^dz] X^I^T^I^ = 8L{B)-'\'/'T'/\ (35) 

\J\z\<l J 



'\z\<l 

Combining the estimates for Ai and A2 gives (|28|) . 

2.2. LEMMA (Delta approximation, part 2). Let ^ e H^(M.'^) and put 



□ 



A 



Then 



z,z)\'dz- / VbA^- dzdz' 



(36) 




<L(fi)-i[Ar-i + 8Ai/^T=^/V/2]. 



(37) 



Proof. Put X{z) = J \'^{z,z')\'^dz', T{z) = J {d.^'^iz, z')\^dz' . By (|8]) we have 
\^{z,z)\'- [ VBA^-^')\'^iz,z')\'dz' 



< L{B)-'[X{z)r-^ + 8X{zy/^T{zf^r^/% (38) 

Integration over z, using the Holder inequality to estimate J X{zy^'^T(zY^^dz, 
gives (0). □ 
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3 Upper Bound 

Let ip G iS(]R^^) be a smooth and rapidly decreasing wave function in the lowest 
Landau level at field strength 1, and let G C^(R^). If ip and 4> are normalized, 
i.e., IV-P = J^N |0P = 1, then 

^b{x^,z) = {BL{B)f/^^{B^/\-^)(f){L{B)z) (39) 

is a normalized wave function in the lowest Landau band at field strength B. 
Moreover, using ( ]I5| ) and (|1^) we have 

E{B,Z,N) < {^b,Hz,n^b) 



where ^iy^) is given by ([18|). Since L{B)'^ / {In B)^ — > 1 as i? ^ oo and ip is 
normalized, one has for the upper bound in Theorem |1.1| only to check that 

and 

6{zi-z,mz)\''d''z 



as 5 — > oo. But this is taken care of by Lemmas 2J. and |2.2| . (That VB,r(^) is 
not defined for r = is of no consequence here, because the error terms in (^) 
and ( pTf ) are integrable all the way to r = 0.) We therefore have 



3.1. PROPOSITION (Upper bound). 

Remark. It is clear that our upper bound holds for fermions, although e{Z, N) 
is the bosonic ground state energy of (0). In fact, in the ansatz (^) above we may 
choose to be antisymmetric and (p to be symmetric; then ^E's is antisymmetric. 
Note also that for the Hamiltonian (|^) the bosonic ground state energy is the 
same as its ground state energy without symmetry restriction. 
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4 Superharmonicity 

In this section we take a closer look at the dependence of the ground state energy 
Ez,n{x'^) of the Hamiltonian (|Hp on the parameter x-^. We start with a simple 
estimate: 

4.1. LEMMA (Simple bounds). The function x"*- ^ Ez^n^x^) satisfies the 
bounds 

N 

- E^' + [sinh-i((Z|x,^|)-i)]') < EzA^^) < (41) 

i=l 

on the set 

A = {x^ e M^^ : xf ^ 0, for all i = 1, . . . , N}. (42) 

Proof. The non-positivity of E is straightforward from the definition by an ap- 
propriate choice of Note that this also holds in the case where some of the 
xf- variables coincide. The lower bound on Ez,n{x'^) follows from Lemma 2.1 in 



LSYall together with the operator inequahty 

Z 



N / 



which is obtained by ignoring the positive two-body interactions. □ 

Next we turn to the superharmonicity properties of Ez,n{x'^)- We shall need 
the following general result. 

4.2. LEMMA (Inherited superharmonicity). Let U be an open set in R.^ 
and assume that / : [/ x M ^ (— oo, oo] is a superharmonic function with the 
property that 

b = minjliminf f{x, t), liminf f{x, t)} 

is independent of x for all x E U. Then g{x) = mftf{x,t) is a superharmonic 
function on U . 

Proof. We shall prove this by showing that /\g < as a distribution. We shall 
use that / is a lower semicontinuous function satisfying the mean value inequality 



f{y,s)dyds < /(x, t)Q+ir'^+\ 

\{x,t)-{y,s)\<r 



for all {x,t) G f/ X R if r > is small enough, where Cd+i is the volume of the 
unit ball in R°'+^ 
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For X G f/ it follows from the lower semicontinuity of / that we have either 
g{x) = b OT there exists t G M such that g{x) = f{x,t). In the first case we 
obviously have 

Cd+ir'^+'gix) > 2 [ g(y)^r^-(x~y)^dy (43) 

J \x—y\ <r 

since g{y) <b for all y. If g{x) < b we also conclude the above inequality since 
^7(a:)crf+i/+^ = f{x,t)cd+y+'> [ f{y,s)dyds 

J\{x,t)-{y,s)\<r 

> / g{y)dyds = 2 g{y)Vr^ - {x - yfdy. 

J \{x,t)~{y,s)\<r \x—y\<r 

Note now that for any G C^{U) we have for any x G t/ that 

lijn^-(d+3) f [(/,(^) - 0(x)] -{x- yfdy = CA0(x) 

'■^O J\x-y\<r 

for some constant C > and in fact this limit holds in the topology of C^{U). 
Thus if (/) > we have 



g{x)/\(t){x)dx = C-^ limr-('^+3) / ^r^ - {x - yfdydx < 

^-^0 J\x~y\<r 

by the inequality (^). Hence A^f < 0. □ 
4.3. THEOREM (Superharmonicity of the energy). On the set A defined 



in ^4M ) ihe function x t— > Ez,n{,x ) is superharmonic in each of the variables 



xf, i = 1, . . . , N independently. 

Proof. We follow closely the proof of Prop. 2.3 in ||LSYa]|, which stated the su- 



perharmonicity of the ground state energy of a one-body operator which can be 
considered as a mean field approximation of Hz,n{x'^)- 

It is clearly enough to prove that Ez^n{x'^) is superharmonic in x^ (on the 
region Xi ^ 0) for x^, . . . , xj^ fixed. We shall prove this by showing that xj- i— *• 
Ez,n{x'^) satisfies the mean value inequality around any given point xj^Q. Let 
Xq = {xj-Q,X2,---,xj^). Choose a sequence of normalized functions G 
C^(M^) such that ifz,7v(x^)^n) ^ ^z,iv(x^) as n -> oo. 

For ti? G M denote by '^l^'' the function 

We clearly have 

inf Hz.N{d)'^^n^) ^ EzAilo) as n ^ oo. 
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If Xi is close to we shall use "^^"^ as a trial function for H[x_-^). We then 
obtain 

EzA^^) < liminf inf (^(^^ ^z,7v(x^X"^ 

Hence 

i^z,iv(x^)-i?z,iv(xo^) < liminf [inf (vl/H, i7z,iv(x^X'^)) - inf (vl/(r), if^,^(xo^)vI/M) 

(44) 

The potential appearing in Hz,N{iI±~), i-e., 

N ^ N 



is a superharmonic function of {zi, Xi) G M'^ \ {0}. Writing 

we see that (^'i'"^ PV^_Ar,j.x\E' is superharmonic in {w,Xi) away from the line 

Xi = 0. Since {'^^\ dl.^^n^) is independent of w and a;]^ for alH = 1, . . . , we 

have that {^^n\Hz,N{d!±~)'^n''') is superharmonic in [w,Xi) away from the line 
xi = Q. 

Moreover, we also have that the two limits 

liminf(^("'),i7z,^(x^X"')) 

are independent of x^. This is true simply because the contribution from the 
terms in the Hamiltonian depending on x^ tend to zero as — > ±cxd. We may 
therefore apply the above lemma to the function f{w, x^) = {'^^\ ifz,Ar(ai"'")^n"'''). 
We conclude that the function 

xi^ mU¥:'\HzA!!i-')^'n ^) 



well 



is superharmonic for 7^ 0. Moreover by the inequality (41) this function is 
bounded below if \xi\ is bounded away from 0. 

Now using Fatou's Lemma we see from (|^) that the average of Ez,n{x'^) — 
En{x^) over the set {xj- : \xi — x^qI < ^} is non-positive for all r > small 
enough. □ 



5 Lower Bound 

The first lemma in this section concerns the ground state energy ef ^(y-*") of 
^z.Niy^) does not use superharmonicity. 
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5.1. LEMMA (Lower bound on ef ^(y"*")). Let K, be a compact subset of the 
set A given in ^4^). Then 



liminf inf e|jv(2/ ) > e{Z,N). 



(45) 



Proof. To avoid problems at points y with — y 

iz.. 



■ for some i,j, we replace 

the repulsive potential VQ^\y±_y±\{zi—Zj) by the smaller potential VQ^\y±_y±\_^_i{zi — 

Zj). We denote the corresponding Hamiltonian by ^iy^) and its ground state 
energy by ef ^(t/-*-). It is obvious that j^f{y^) > ef ^(y-*"), so a lower bound on 
^z,N{y^) gives a lower bound on ef^^d/-*-). 

Let be a normalized, symmetric wavefunction in C(j"(M^). Since {"^j hz,N^) > 
e(Z, N) we have to estimate the matrix elements of the difference ^{y-^) — hz,N- 



Using Lemmas and |2]^, together with the Holder inequality for the integration 
over Z2, . . . ,zn and z^, . . . , zn respectively, we obtain 



hl^^) - hz^N"^) < L{B)-\ZN + N{N - 1)) 



X 



^mL + 8T^/"(2w + 1 



>l/2 



(46) 



where rmin and rmax are respectively the minimum and the maximum value of 



N, with y G /C, and 



= N J \dz^{z, Z2, . . . , ZN)\'^dzdz2 - ■ ■ dzN 



(47) 



is the kinetic energy of \E'. Now if \E'"", , n 



1,2,... is a minimizing sequence 

of normalized wave functions for j^{y^), then we may assume that the corre- 
sponding kinetic energy is uniformly bounded in n, B and y-^ G /C. In fact, we 
may assume that (\E'^x „, Ar(l/"'")^^x „) is a bounded sequence. If we use the 
bound from Lemma 2.1 in [LSYa], we obtain 

(^fx„,/^f,^(y^)%^x J > y-^ (^)' (l + [sinh-{(2Z)-5V^}]^) , (48) 

where we have saved half of the the kinetic energy T„ of . For large B, 

y_ i'^ 

L{B)-^ [sinh-^{(2Z)-i5i/2}] is bounded and hence we see that T„ is bounded. 
The error term (^) with \I' = thus tends to zero as 5 ^ cxd, uniformly in 



n, and the lemma is established. 



□ 



5.2. LEMMA (Uniform bounds on Ez^n{x^))- Let e > 0. Consider the set 



C 



B,e 



{x^ : £5-1/2 < 



for alli = 1,...,A^}. 



(49) 
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Then 



limmf(lnS)-^inf{Ez,7v(^^) : e C^'"} > e{Z,N). (50) 



where e{Z, N) as before denotes the 1- dimensional delta function atom energy. 

Proof. Define the sets 

= {x-^ : eB-^/^ < \x:t \ < n, for alH = 1, . . . , N}. 

Since C^'^ is compact and Ez^n is lower semicontinuous (being superharmonic, 
in fact, superliarmonic in eacli variable) we may find xj^ G C^'^ such that 

Clearly, 

lim EzMl^in) inf{^z.^(x^) : G C^'^}. 

n— >oo 

By the superharmonicity of Ez,n{x'^) in each variable xj- we know that each 
coordinate of the point satisfies either = e_B~-^/^ or = n. 

Moreover, since Ez,n{x'^) is invariant under permutations of the coordinates of 
X"*- we may assume that |a;^„| < < . . . < \xj^n\ ^- By possibly going 

to a subsequence we may assume that there exists an integer < K < N such 
that for n large enough 



\x^ I 



£5-1/^ for2 = 
n, for i > K 



Moreover, we may assume that x/-^ converges as n ^ oo for i = 1, . . . , i^. 

Since we may ignore the variables x^„, i = K+1, . . . , N, which tend to infinity 
we have 

lim Ez,n{x^)/Ez,k{xI^, xj^ ) = 1 

Since Ez^xix-^) is lower semicontinuous we conclude that there exists a point 
e M^i^ with |x,^^| = for alH = 1, . . . , such that 

inf{Ez,7v(x^) : x-^eC''''} = Ez,K{^too,---,xioo)- 



By Lemma |5.1| we have that 
liminfinf{L(5)-2i?^,,,^(fi-V2^^^,... ,5-1/2^^.) ; \y^\ = e, for alH} > e(Z, iT). 

Since K < N and hence e{Z, K) > e{Z, N) we have proved the lemma. □ 
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5.3. LEMMA (Wave functions in the lowest Landau band). If-^eH% 
belongs to the lowest Landau hand at field strength B, then J {"^{x-^, z)\'^dz is a 
bounded function of x'^ (possibly after a modification on a null set) and for all 
l<n< N 



sup 




\'^{x'^,z)\^dzdx^_^_-^ ■ ■ ■ dxjf 



TDn 



Proof. The projector 11^ on the lowest Landau band is the A^-th tensorial power 
of the projector 11° that operates on L^(M^; C^) and is given by the integral kernel 

nO(x, x') = Ul{x^, x'^)6{z - z')P^, (52) 

where 

n'[(x^,a;'^) = ^exp \ Ux^ x x'^) ■ B - - x'^)^b} (53) 

and is the the projector on vectors in with spin component —1/2. The 
kernel n*[(x^, x'"*") is a continuous function with 



{x^,y^) (54) 



and 

U'{x^,x^) = ^ (55) 

for all a;-*". A wave function in the lowest Landau band has the representation 
= n^\E'. After writing 11^ as an integral operator ( pT| ) follows from the Cauchy- 
Schwarz inequality, using (|5^) and (^5[). □ 



5.4. PROPOSITION (Lower bound). 

liminf^^Mli^>e(Z,A). (56) 

Proof. For fixed B let \E' be a normalized wave function in the lowest Landau 
band. By (|T5D we have 

{^,Hz,n'^) > J Ez,n{x^) (^j \'^{x^,z)\^d^ dx^. (57) 



We split the integral over x_ into an integral over C (defined in (|49|) ) and its 
complement in M^^. By Lemma |5.2| we have only to consider the latter. Using 



the estimate (^Tj) the task is to bound terms of the form 

;i + [sinh-^(Z|x^^|-^)]2) (^j \^{x^,z)\^d^ dx^ (58) 



xi\<eB-^/^ 
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from above. If i = j we carry out the integration over all with k ^ i 
and use Lemma |5]^ for the remaining variable x^. For small r, | sinh~^r^^| < 



(const.) I lnr| and the term can be estimated by 

(const.)/ {\n\x^\)^Bdx^ <{const.)e^{\nBf. (59) 

For i ^ j we split the integration over xf into two parts, namely < B^^^"^ and 
\x-j \ > 5^^/^. For the first part we obtain the following bound, after transforming 
variables and using Lemma 5.3, this time for n = 2, 

(const. ) /" {\YiB'^/'^\yf\fdyidyf<{consi.)e^{\nBf. (60) 

J\yi\<Uyf\<'^ 

For the integral over \x-^\ > B^^^"^ we estimate | sinh~"^(Z|xj"|~"^)p by its maxi- 
mum value, < (const.) (In 5)^ and obtain for this part of the integral the upper 
bound 



ix , z)\'^dz I dx 



xf\<eB- 



(1 + (const.) (In 5)2) (^j \m{ 
< [ (1 + (const. ){\nBf)Bdxi < e^{l + c{\nBY), (61) 



where we have used Lemma ^]3| again. We see that (55) is bounded above by 
(const.) (e In 5)2, for B large enough. Since e > is arbitrary this completes the 
proof. □ 



6 The one-dimensional delta-function model 

We now want to study the the delta-function Hamiltonian (|^, in particular its 
mean field limit, N ^ oo, Z ^ oo, with A = N/Z fixed. 

For this it is convenient to make a scale transformation z z/Z, which 
implies a unitary equivalence 

hz,N = Z'hz,N (62) 

with 



hz,N = J2 (P' - '^(^^)) + I E ^('^ - '^)- (63) 

i=l i<j 

We denote its ground state energy (again in the sense of quadratic forms) 
by e{Z,N). The formal mean field theory of this system is identical to the so 
called hyper-strong theory discussed in ||LSYa|| , Section 3. The energy of a (one 



dimensional) electron density Zp in this theory is [p] with 



BSY 19/11/99 



16 



^^'[p] = / t: VpM - P(0) + / p'dz (64) 



The infimum over densities with fixed normahzation j p = X leads to the hyper- 
strong energy E^^{X) given by the right side of (^. 

We shall now establish this mean field limit rigorously and prove Theorem 

1.2. 

6.1 A comparison model 

An upper bound to the Hamiltonian (|^) can be obtained from another model 
whose ground state can be computed explicitly. The corresponding Hamiltonian 
is completely symmetric with regard to each single refiection Zi — > —Zi, and the 
electronic repulsions are equally distributed between the sites Zi and —zf. 

^ 1 
hz,N = J2iP^- '^(^^)) + ^ 5Z - ^^•) + '^(^^ + ^^•)] • (65) 

i=l i<j 

Its ground state energy is denoted by e{Z, N). 

The replacement of 1/Z by 1/{2Z) is important, because it compensates to a 
certain extent the doubling of the interaction sites. In particular it leads to the 
same formal mean field theory as ( |63l) for symmetric electron densities p. This 
observation will be substantiated by the mathematical treatment in the sequel. 

The model ( |65|) was used in |[WS|| for = 2 as a starting point for a perturba- 



tional calculation. It was also considered in |[Ro|| (for = 2) as an upper bound 
to the model (|^), but with the coupling 1/Z instead of 1/{2Z). The present 
considerations and extensions to A^ > 3 appear to be new. 

The ground state wave function ip , Hit exists, is completely symmetric under 
permutations of {zi...zi^} and refiections Zi —Zi. Such a highly symmetric 
function ip is determined by its restriction to the cone 

M = {z: < zi < Z2 < . . . < Zn} (66) 

In A4 we make the ansatz 

N 

^{z,...zr,)=cl[e-^^'% (67) 

i=l 

with c a normalization constant and let hz,N act on the symmetrically extended 
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The delta-function interactions dictate the jumps in the partial logarithmic 
derivatives of ip at the boundary of Ai and we find 

_ 1 _ 1 

2 AZ 

which implies 

1 n — 1 

^-=2-^- ^^^^ 
The function (|67D is square integrable if and only if all are strictly positive, 
which is equivalent to 

N<2Z + 1. (69) 

The corresponding eigenfunction ip, is everywhere positive, and it is easy to 
see that it is, indeed, a ground state for (|65D: Define the operators An on L^(M^) 
by 

A. = 9,„-9,„(ln^), (70) 

with obvious domains of definition. Denoting by e('?/') the eigenvalue of hz,N 
corresponding to we can write the quadratic form hz^N as 

N 

hz,N = J2'^nAn + e{^)- (71) 



i=l 



The equation 



N 



{i^,hz,Ni^) = II ^-i^ II' II ^ II' ^ ^(^) II ^ II'' (72) 



i=l 



which holds for each ip in the form domain of hz,N, shows that e{Z, N) = e{ip). 

If N > 2Z + 1, the simple inequality e{Z,N) < e{Z,N — 1) is sufficient 
for our purposes. To prove it, one may use trial- wave-functions of the form 
ip{zi . . . ZN-i)£f{£^ZN) whith a smooth if, and take e to zero. This inequality for 
the energies can be iterated to 

e(Z,iV) <e(Z,iVo), (73) 



where is the largest integer satisfying 

For < 2Z + 1 the ground state energy is the eigenvalue corresponding to 



N 



n=l 
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li N > No we may use (73), i.e., e{Z, N) is bounded from above by (|73|) with A 
replaced by Ao = Nq/Z. Dividing by Z and keeping A fixed, the leading term for 
— > cxD is identical to E^^{X). By the next proposition this is sufficient for the 
upper bound in Theorem 1.2. But one can in fact show that e{Z, No) is equal to 
e{Z, N) for iV > 2Z + 1 and not only an upper bound to it. Hence No is equal to 
Nc, the maximal number of electrons that can be bound in the model (|65D. We 
give the proof of this result in the appendix. 

6.1. PROPOSITION (The comparison model gives upper bounds). The 

ground state energy of the symmetrized model hz,N is an upper bound to the 
ground state energy ofhz,N'- 

e{Z,N) <e{Z,N). (75) 
This inequality is strict, if N >2 and N < 2Z + 1. 

Proof. The Hamiltonian Hz^n is the symmetrization of hz^N with respect to the 
group 71 with 2^ elements, generated by the reflections Zi —Zi, i = 1, . . . ,N. 
For i? G 7^ let Ur denote the corresponding unitary operators on L^(]R^) . Then 

^ '^{URi^,hz,NURij) = {lpJlz,Nij) 
ReTZ 

for any so 

e{Z,N) <e{Z,N). 

If N < 2Z + 1 we may take the square integrable ground state wave function 
of hz,N, given by (|67|) , as a test state for hz,N- It satisfies UriP = for all R, so 

{^,hz,N^) = e{Z,N). 

But ifj is not an eigenfunction of hz^N ii N > 2, so e{Z, N) is strictly below 
e{Z,N). ' □ 

Combining the last proposition with Eq. (|74|) , recalling that e{Z, N) = Z^e{Z, N), 
we obtain 

6.2. PROPOSITION (Upper bound in the mean field limit). IfN,Z ^ 
oo with A = N/Z fixed, then 

limsupe(Z, A^)/Z=^ < ^"^(A) (76) 

where E^^{X) is given by (j^. 
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6.2 Lower bounds to the delta-function Hamiltonian 

An elegant way to obtain lower bounds for Hamiltonians with repulsive pair 
interactions is the use of positive definite functions. This was probably done for 
the first time in [pdL'l]| . In this method, the positive definite functions have to be 
finite at the origin, however, and hence it is impossible to bound the 5-function 
interaction in this way without additional help. Our way out is to borrow a bit 
of kinetic energy (this was also done in Theorem 7.1 in ||LSYa|| ). So we search for 
operator inequalities 

ap^ + -6{z)>wz,a,b{z) (77) 

with appropriate functions Wz,a,b{z), depending on a parameter b in addition to 
a and Z to allow convergence to a delta function. 



or 



6.3. LEMMA (An operator inequality). The inequality ( f7?| j holds ft 

Proof. With the simple reformulation to 

ap^ + ^6{z)-wz,a,b{z)>0 (79) 

we are on well known territory: The Hamiltonian on the left side shall have no 
negative eigenvalue. By the scale transformation z Zaz, this inequality is 
transformed to 

/ + 6{z) - Wb{z) > 0, Wb{z) =Z^a wzM^az). (80) 
This inequality will hold for 

if it is true for the larger potential 

_ A2 -b\z\ 

Wdz) 



(26 + l)(l-e-^l^l 



because the Hamiltonian in ([80| ) is bounded from below by 

p^ + 6{z)-W,{z). (82) 

This Hamiltonian has 
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as a positive symmetric solution to the Schrodinger equation - as a differential 
equation - with zero energy. 

Now, if ( p2D would have a square integrable ground state wave function g{z), 
this wave function would also be symmetric under reflection z — —z, and the 
delta-function would dictate the same value for g\z)/g{z) as it does for f\z) / f{z) 
at 2 = 0+. So the question of the existence of g{z) can be dealt with by the 
methods which are used for proving Sturm's comparison theorem: We assume 
that g{z) exists, with negative energy E. The Wronskian W{z) := f'{z)g{z) — 
g'{z)f{z) is zero ai z = 0+. Its derivative is determined as W'{z) = Ef{z)g{z). 
If g{z) is chosen positive, then W'{z) is negative, which implies that W{z) is 
negative for z > 0, and g'{z)/g{z) > f'{z)/f{z). This inequality can be integrated 
to give g{z)/g{0) > f{z)/f{0), a contradiction to the assumption of the square- 
integrability of g(z). 

Therefore we know that the Hamiltonian (^) has no negative eigenvalue. 
And so the operator inequality holds. □ 

The {VFf)(-2)} and hence {Zwz,a,b{z)} are 5-sequences in the limit b ^ oo. All 
these functions are positive definite, and finite at the origin: 

wz,a,bm < (84) 

With this tool we can now deduce the lower bound for the many body Hamil- 
tonian: 

6.4. PROPOSITION (Lower bound in the mean field limit). IfN,Z ^ 
oo with A = N/Z fixed, then 

\immfe{Z,N)/Z^>E^^{X). (85) 

Proof. We use the operator inequality ( 177|) with wz{z) := wz,a,b{z), (a and b will 
finally be chosen as appropriate powers of A^) to bound hz,N from below. For 
each 6{zi — zj) we use it twice; one time with apf, and a second time with ap'j. 
Then we add these inequalities and divide by two: 



N 

E 

1=1 

AT 



I — a- 



N -1 



I — ^ + -6{zi-Zj^ 



(86) 



i=l 



i<j 



At this point the positive definiteness of wz{z) becomes essential. It implies, 
that for any real valued integrable function criz): 



> 
(87) 
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Expanding this expression and integrating the delta-functions we get 



Ar2 



wz{zi - z)a{z)dz - —wz{^) 



~Y I I ^^^)^z{z - y)a{y)dzdy. 



Combining this with ( pUf ) gives 

hz,N >^hi{Z,N,a) — / / a{z)wz{z - y)a{y)dzdy (89) 

i=i ^ ^ 

with the one-particle operators 

hi{Z, N, a)=(^l- a^^^) Pi - S{z,) + N{a * wz){zi) - ^wz{0). (90) 

The parameters are now chosen as 

a = N-^-', b = N', with 0<£< 1/2. 

The fraction of kinetic energy per particle that we borrowed in (|86D then decreases 
as A^~^, and the functions wz{z) become 

^ Z2 (2iV- + l) • ^ ^ 

In the mean field limit N, Z ^ oo with N/Z = A > fixed the sequence Zwz{z) is 
a (5-sequence, and wz{0) ~ XN'^'^ /Z 0. If (t{z) is smooth with |a"'(2;)| < 7, then 
\N{a*wz){z)—cr{z)\ < 2'-)\^N~^ . The one particle Hamiltonians h{Z, N, a), with 
smooth cr{z), converge as quadratic forms pointwise (i.e., for each test function) 
to 

hxa = P^-5iz) + Xa{z). (92) 

Moreover 

h{Z, N, a) > hxa - {N-'/2)p^ - 26X^N~' - {X^ /2)N^^~\ (93) 

Since the ground state energies of operators of the type ap^ + V are concave 
functions of a and hence continuous in a, the ground state energies of the right 
side of ( p3D converge in the limit 00. 

The ground state energy of h\„ is a concave functional e[Acr], and the lower 
bound (P^, when divided by the number of electrons A^, gives 

limmf ie(Z, N) > e[Xa] j a\z)dz =: IM- (94) 

Inserting the mean field density p for Act (i.e., the minimizer of (64) which satisfies 
Eq. (3.8) of ||LSYa|| ) gives the mean field energy, divided by A, as a lower bound 
to the limit of the energy per electron. □ 
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We remark that searching for the supremum of I\ [a] in (Q) also leads to the 
mean field equation of |[LSYa|| : Assuming e{\(j) = {ip, hx^ip) with a normalized ifj 
the variational condition on a{z) for maximizing T\[ct] is 



(95) 



Inserting this into the Schrodinger equation /iao-"^ = /WaV' gives 

- i)'\z) - 5(^)V^(0) + \i)\z) = -fix^iz), (96) 



i.e., Equation (3.8) in [|LSYa| . 

Finally we remark that the energy per electron, e'{Z,N)/N, approaches the 
mean field limit monotonously. There is also a subadditivity property, which in 
the limit becomes concavity of E^^{X)/X. These properties of the approach to a 
mean field hold in some other cases too, as will be shown elsewhere P99|| . 



7 Conclusions 

We have shown that the energy of an atom in a strong magnetic field B ap- 
proaches, after division by (In 5)^, the energy of a many body Hamiltonian with 
delta interactions in one dimension as S — >■ cx). This delta function model is not 
explicitly solvable, but an upper bound to the energy can be given in terms of 
another model with the same mean field limit and where we can explicitly calcu- 
late the ground state energy. In the latter model an atom with nuclear charge Z 
can bind up to 2Z electrons. Whether this represents the true state of affairs for 
the atomic Hamiltonian in the B oo limit is an open problem. 
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Appendix 

We prove here that the ground state energy e{Z, N) of the Hamiltonian (|65|) is 
independent of if > 2Z + 1. 

PROPOSITION (Maximal negative ionization for the comparison model). 

IfN>2Z + l, then 

eiZ,N)=eiZ,N,) (97) 
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where No is the largest integer strictly smaller than 2Z + 1. Moreover, there is 
then no L?' -function with e{Z, N) as an eigenvalue. 

Proof. In the cone A4 defined by (pBI) we consider the wave function 



tP{z,,...,z^) = l[e-^''^ II (l-n^z,), (98) 

i=l j=No+l 

with Kn defined by ([681). Since kj < for j > + 1, the function ip is strictly 
positive. We extend ip symmetrically from A4 to all of as a continuous 
function. 

The jumps in the logarithmic derivatives of tp at the boundary of Ai are not 
of the right size required for an eigenfunction of hz,N- But tp is an eigenfunction 
of a slightly different operator: 

hz,Nip = e{Z,N)tP (99) 

with 

No 

eiZ,N) = -Y,^l = e{Z,No) (100) 
1=1 

and 

^ 1 

hz,N = ^{p1- H^i)) + -^^lijizi, ...,Zn) [S{zi - Zj) + S{zi + Zj)] . (101) 

1=1 i<j 

with certain functions 7jj. It is sufficient to specify 7,^4+1 (zi, za?) on the 
boundary of M (other cases follow by permutation and/or reflection of the vari- 
ables) and one finds for < zi < 2:2 < . . . < zn'- 

( 1 a l<i<No-l 

7i,i+i = I AZ {kno + l«^7Vo+i|(l + \i^No+i\zNoy^) a i = No . (102) 

( {1 + \tii\zi)-\l + \Ki+i\zi+i)'^ if No + l<i<N 

Since 7j,j+i < 1 for all i one has 

hz,N < hz,N. (103) 
Since ip is stricly positive we can in the same way as in (^) write 

N 

hz,N = Y.AlAn + e{Z,N) (104) 

n 

with 

in = a.„-9,„(ln^), (105) 
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and conclude that e{Z, N) = e{Z, No) is, indeed, the ground state energy of hz,N- 
Hence, e{Z, N) = e{Z, N^) ioi N > 2Z + I. 

To see that there are no bound states at the bottom of the spectrum of hz,N 
assume ip is an eigenfunction to eigenvalue e{Z, N), so that 

{^p,hz,N^) = eiZ,N)\\^pf. (106) 
By ( |103| ) and the equality of the ground state energies this implies 

{tP,hz,N^) = e{Z,N)\\^\\\ (107) 
which, because of (|104|) , is equivalent to the set of differential equations 

i^V^ = 0. (108) 
These equations have no other solutions than op, and tp is not an function. □ 
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